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NOTATION 


Area 

Magnetic  flux  density  (or  magnetic  induction)  vector 
Velocity  of  light  in  free  space 
Electric  displacement  vector 

Electric  field  intensity  vector  (in  original  medium) 

Electric  field  intensity  vector  in  model 

Ratio  of  model  characteristic  electric  field  intensities  to  those  in 
original  medium 

Frequency  (and,  including  and  following  Equation  (10c),  ratio  of 
model  to  original  frequencies) 

Magnetic  field  intensity  vector  (in  original  medium) 

Magnetic  field  intensity  vector  in  model 

Tangential  component  of  magnetic  field  intensity  vector 

Ratio  of  model  characteristic  magnetic  field  intensities  to  those  in 
original  medium 

/“ 

Free  electric  current  density  vector 
Propagation  constant  of  light  when  o  =  0 
Characteristic  length  (in  original  system) 

Characteristic  length  in  model 

Ratio  of  model  characteristic  lengths  to  those  in  original  system 

Induced  magnetization  vector 

Induced  electric  polarization  vector 

Time-average  electromagnetic  power  absorbed 

Mechanical  scale  factor  (reciprocal  of  1) 

Poynting  Vector  (S  =  E  «  H) 

Ratio  of  model  characteristic  times  to  those  in  original  system 


Characteristic  time  (in  original  system) 
Characteristic  time  in  model 


Velocity  of  light  in  a  material  medium  (v  =  1 //em) 

Rectangular  Cartesian  coordinates  of  a  point  in  space 
Complex  characteristic  (or  wave)  impedance  of  a  medium 
Real  part  of  the  characteristic  impedance  (resistance) 

Imaginary  part  of  the  characteristic  impedance  (reactance) 
Characteristic  impedance  of  a  medium  when  o  =  0 

Real  part  of  the  complex  propagation  constant  or  attenuation  constant 
Imaginary  part  of  the  complex  propagation  constant  (=  k  when  o  =  0) 
Complex  propagation  constant 
Skin  depth 

Partial  derivative  operator 

Electric  permittivity  (dielectric  constant)  in  original  medium  (and, 
following  Equations  (13a)i  ratio  of  model  to  original  permittivity) 

Electric  permittivity  in  model 

Electric  permittivity  of  free  space 

Ratio  of  model  characteristic  impedances  to  those  in  original  medium 
Electric  loss  tangent 

Electric  permittivity  relative  to  free  space 
Magnetic  permeability  relative  to  free  space 
Wavelength 


Magnetic  permeability 
(13a),  ratio  of  model 

Magnetic  permeability 

Magnetic  permeability 

Pi  (3.141592654) 

Fri  e  cha'’ge  density 


in  original  medium  (and,  following  Equations 


in  model 
of  free  space 


Accesion  For  I  "  ] 

ntis  cra&i 

DTIC  TAB 

Unannounced 

Justification 

^ — 

□ 

□ 

.  . . . 

By . 

Di:4  ibution/ 

- - 

Availability  Codes 


0 


Electric  conductivity  in  original  medium  (and,  following  Equations 
(13a),  ratio  of  model  to  original  conductivity) 


Electric  conductivity  in  model 
Electric  susceptibility 
Magnetic  susceptibility 
Angular  frequency  (=  2Trf) 
Divergence  operator  (div) 

Curl  operator 

Laplacian  operator  (div(grad)) 
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ABSTRACT 


This  report  presents  a  complete  and  rigorous  treatment 
of  the  predictions  of  electrodynamic  similitude  relating 
scattering  from  scaled  targets  to  scattering  from  full-size 
targets.  The  usual  assumption  is  made  that  the  targets 
treated  are  those  whose  conductivity,  permittivity,  and 
permeability  do  not  change  with  frequency.  Three  general 
nonlinear  modeling  equations  in  six  variables  are  derived, 
and  a  complete  set  of  solutions  is  presented.  For  the  first 
time,  particular  emphasis  is  given  to  the  effects  of 
differences  between  geometric  and  complete  scaling  on  the 
electromagnetic  fields  and  on  the  radar  cross  section,  and 
effects  of  approximations  to  complete  scaling  are  evaluated. 
Conditions  are  obtained  on  properties  of  materials  required 
for  models  made  from  these  materials  to  accurately  simulate 
systems.  Absorption  and  energy  balance  are  also  treated, 
and  the  influence  of  finite  conductivity  on  surface  currents 
is  shown. 

The  possibility  of  using  scaled-down  models  of  real 
targets  is  usually  claimed  to  be  based  on  the  linearity  of 
Maxwell's  equations.  Therefore,  we  first  examined  the 
physics  involved  in  linearity,  its  relation  to  the  Maxwell 
equations  in  macroscopic  media,  and  its  impact  on  physical 
scale  modeling.  The  behavior  of  electromagnetic  fields  was 
also  examined  as  a  function  of  frequency.  Upper  bounds  on 
the  electric  conductivity  and  limitations  on  the  electric 
and  magnetic  polarizability  of  realistic  materials  may  have 
profound  implications  for  model  measurements.  All  the 
energy  density  in  the  reflected  wave  may  be  lost  at  higher 
frequencies  by  incomplete  scaling  of  even  apparently 
nonabsorptive  targets.  Thus,  the  ratio  of  measured  radar 
cross  sections  of  geometrically  scaled  targets  to  those 
predicted  for  a  completely  scaled  model  will  eventually 
approach  zero. 


ADMINISTRATIVE  INFORMATION 

This  project  was  supported  by  the  DTNSRDC  Independent  Exploratory  Development 
Program,  sponsored  by  the  Office  of  Chief  of  Naval  Research,  Director  of  Navy 
Laboratories,  OCNR  300  and  administered  by  the  Research  Coordinator,  DTNSRDC  012.3 
under  Program  Element  62766N,  Task  Area  ZF-66-^1 2-001  under  DTNSRDC  Work  Unit 
1-2930-001 . 


INTRODUCTION 

The  recent  development  of  millimeter  and  submillimeter  and  of  infrared  and 


visible  laser  dinar)  rad  ir  .syst-^ns  for  detection  of  military  vehicles  [1]*  has 
stimulated  oonsideraMe  i  ntf-r-*^;'.’  in  scal“d-down  models  of  radar  targets.  The 
distribution  of  radiation  S'lfte-t^d  from  a  full-size  target  is  expected  to  be 
predictable  from  neasurer''nts  on  a  faithful  scale  model. 

Models  of  electromagnet i^  systems  have  been  used  for  antenna  studies  [2]  since 
the  early  1900's.  However,  the  earliest  known  work  on  model  measurements  of  radar 
cross  sections  (RCS)  was  at  Massachusetts  Institute  of  Technology  Radiation 
Laboratory  in  19^2  [3].  This  application  received  considerable  impetus  from  the 
suggestions  and  work  of  Sinclair  and  his  colleagues  during  World  War  II  [il,  5],  and 
has  been  extensively  pursued  since  then  [6~12].  Despite  some  problems  in  preparing 
the  models  themselves  and  in  making  measurements,  scale  models  provide  an  important 
means  both  of  determining  the  RCS  of  targets  for  which  no  other  method  is  practical 
and  of  providing  independent  verification  of  values  found  in  other  ways. 

Scale  models  of  physical  systems  all  depend  on  principles  of  similitude  which 
allow  relationships  to  be  established  between  model  and  full-scale  behavior. 
Equations  for  scaling  electromagnetic  models  have  previously  been  written  by 
Stratton  [131.  Brown  and  King  [19,  15],  and  Sinclair  [5],  However,  practical 
difficulties  with  scaling  properties  of  realistic  materials  now  exist  which,  for 
most  scattering  measurements,  discourage  all  but  simple  geometric  scaling  [1,  9-12, 
l6-i8].  Geometric  scaling  is  a  partial  form  of  scaling,  where  the  scaled  target 
has  the  same  dimensions  in  scaled  wavelengths  as  the  actual  target  has  in  actual 
wavelengths,  but  material  properties  are  ignored. 

Since  a  complete  set  of  solutions  to  the  general  scaling  equations  has  never 
previously  been  published,  little  concern  was  shown  about  differences  between 
g'^ometric  and  complete  scaling  in  radar  applications  [1].  Although  it  has  been 
Known  that  scaled  electric  conductivity  (o)  should  increase  with  frequency  (f), 

*  References  are  listed  on  page  39. 
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differences  between  surface  current  distributions  on  poor  conductors  and  on 


perfectly  conducting  targets  were  believed  to  be  unimportant.  Thus,  infinite  o  is 
usually  assumed,  the  f-dependence  of  material  properties  is  usually  ignored,  and 
little  attention  has  been  given  to  resulting  differences  in  the  electromagnetic 
fields  involved. 

However,  this  report  will  show  this  lack  of  concern  is  not  justified.  Here, 
we  develop  a  complete  and  rigorous  treatment  of  the  relations  between  scattering 
from  scaled  and  full-size  targets  on  the  basis  of  electrodynamic  similitude, 
starting  with  the  usual  assumption  that  not  only  o  but  also  electric  permittivity 
(dielectric  constant)  (e)  and  magnetic  permeability  (u)  do  not  change  with  f. 

Three  general  nonlinear  modeling  equations  in  six  variables  are  derived,  and  a 
complete  set  of  solutions  is  presented.  For  the  first  time,  particular  emphasis  is 
given  to  effects  of  differences  between  geometric  and  complete  scaling  on  the 
electromagnetic  fields  involved  and  on  the  RCS,  and  effects  of  approximations  to 
complete  scaling  are  evaluated.  Conditions  are  obtained  on  properties  of  materials 
required  for  models  made  from  these  materials  to  accurately  simulate  systems. 
<\bsorpticn  and  energy  balance  are  also  treated,  and  the  influence  of  finite  o  on 
surface  currents  is  shown. 

The  possibility  of  using  sealed-down  models  ot  real  targets  is  usually  claimed 
to  be  based  on  linearity  of  Maxwell's  equations  [6,  10-12].  Thus,  we  first  examine 
the  physics  involved  in  linearity,  its  relation  to  the  Maxwell  equations  in 
macroscopic  media,  and  its  impact  on  physical  scale  modeling. 

The  behavior  of  electromagnet i:  fields  is  also  examined  as  a  function  of 
frequency.  Unfortunately,  it  is  physically  impossible  to  satisfy  scaling 
conditions  exactly  when  realistic  materials  are  used  and  the  change  in  f  is 
apprc' ' iable ,  because  this  would  requi’^e  that  upper  bounds  on  physical  o  bo  '^xc>^cic,i 
at  higher  f,  and  definite  limitations  also  exist  on  physical  e  anc. 
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This  report  will  show  that  incomplete  scaling  of  apparently  nonabsorptive 


targets  nay  have  profound  implications  for  model  measurements.  Both  electric  and 
magnetic  scattering  are  considered.  The  energy  density  in  the  reflected  wave  will 
be  examined  at  higher  frequencies  for  both  geometrically  and  completely  scaled 
targets,  and  the  frequency  dependence  of  the  RCS  measured  for  geometrically  scaled 
targets  will  be  compared  to  that  predicted  for  a  completely  scaled  model. 


THE  MAXWELL  EQUATIONS  IN  MACROSCOPIC  MATERIALS 

All  scattering  characteristics,  including  interference  and  diffraction, 
polarization,  and  creeping  and  traveling  (surface)  wave  phenomena,  should  be 
correctly  represented  by  rigorous  scale  modeling. 

If  the  model  is  to  be  an  accurate  simulation  of  a  full-scale  system, 
transformations  relating  modf^l  to  full-scale  quantities  will  transform  Maxwell's 
equations  from  model  to  full-scale  systems.  Since  electromagnetic  fields  in  both 
systems  must  satisfy  Maxwell's  equations,  we  can  thereby  determine  how  e,  m,  and  a 
must  transform  for  simulation  to  be  accurate.  However,  because  upper  bounds  on 
electric  conductivity  and  limitations  on  electric  and  magnetic  polarizability 
so-i’'t  im<'s  rend-'r  complete  scaling  impossible,  approximations  to  complete  scaling 
s'-',.-  lid  b-:-  inv-^stigated. 

To  examin--'  the  origin  and  impact  of  these  approximations,  we  must  start  from 
thf-  ’I'ixv-'ll  ^qu-ations  in  macroscopic  materials.  Following  the  rationalized 
M---_er-K i  logr  rn-Second-Amp«^re  (MKSA)  system  of  units  used  by  Stratton  [13],  these 


’’lit  ions  can  be  written 
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V  X  E  = - 

3t 


(la) 


n  -  n  _  /^ 


*4.  (i"  ■•j*  .gg*  x*f  ’  J't*'  ^  ft' 
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V  X  H  =  J  + 


V  •  D  =  p 


where 


E  =  electric  field  intensity  (V/m) 

B  =  magnetic  induction  (T  or  Wb/m*) 

H  =  magnetic  field  intensity  (A/m) 

D  =  electric  displacement  (C/m^) 

J  =  free  electric  current  density  (A/m*) 
p  =  free  charge  density  (C/m*) 

and  the  vector  operators  (V  x  and  V  •)  are  used  to  indicate  how  changes  in  space  of 
the  above  quantities  are  coupled  to  their  changes  in  time  (t). 

E  and  B  are  macroscopic  electric  and  magnetic  field  quantities,  averaged  over 
volumes  which  are  large  compared  to  the  volume  occupied  by  a  single  atom  or 
molecule.  D  and  H  are  corresponding  derived  electric  and  magnetic  fields,  related 
to  E  and  B  through  constitutive  relations 
D  =  d[e,  b] 

(3) 

H  =  h[e,  b] 

Square  brackets  are  used  to  signify  functional  relations  that  are  net  necessarily 
simple,  because  they  may  depend  on  past  history  (hysteresis),  may  reflect  the 
derived  polarization  fields'  lagging  behind  the  changes  in  source  fields  causing 
these  effects,  may  be  nonlinear,  or  some  other  factor. 

The  two  homogeneous  equations.  Equations  (la)  and  (lb),  can  be  solved  formally 
by  expressing  E  and  B  in  terms  of  scalar  and  vector  potentials,  but  the  two 
inhomogeneous  equations.  Equations  (2a)  and  (2b),  cannot  be  solved  until  the 
derived  fields  D  and  H  are  known  in  terms  of  E  and  B.  Bound  charges  and  currents 
enter  Equations  (3)  through  these  derived  fields,  which  represent  macroscopically 


-•.jr.y.-r 


averaged  electric  and  magnetic  dipole,  electric  and  magnetic  quadrupole,  and  higher 


moment  densities  of  material  media  induced  by  applied  fields. 

In  addition,  forces  on  electrons  depend  on  both  magnetic  and  electric  fields, 
so  for  currents  in  conducting  materials  there  is  a  generalized  Ohm's  law, 

J  =  ^[e,  b]  (4) 

In  most  materials,  only  electric  and  magnetic  polarizations  P  and  M  are 
significant.  Then,  electric  quadrupole  and  higher  terms  in  Equations  (3)  are 
completely  negligible,  and  we  may  define 


D  =  GqE  +  P 
B  =  Po(H  +  M) 


(5) 


where 

P  =  induced  electric  polarization  (C/m^) 

M  =  induced  magnetization  (A/m) 

£o  =  permittivity  of  vacuum  =  8.85^  pF/m 
Mo  =  permeability  of  vacuum  =  400it  nH/m 
These  equations  form  the  basis  of  all  classical  electromagnetic  phenomena. 
When  they  are  combined  with  the  Lorentz  force  equation  and  Newton's  laws  of 
classical  mechanics  to  describe  coupling  between  the  fields  and  their  moving 
sources,  they  provide  a  complete  description  of  the  classical  dynamics  of 
interacting  charged  particles  and  electromagnetic  fields  [19]. 


IMPACT  OF  LINEARITY  ON  PHYSICAL  SCALE  MODELING 
The  claim  that  scaled-down  models  can  be  used  to  simulate  real  targets  is 
based  on  linearity  of  Maxwell's  equations  [6,  10-12].  Hence,  our  treatment  does 
not  apply  to  obviously  nonlinear  materials,  such  as  ferromagnetic  materials  and 
ionized  regions  with  magnetic  fields.  However,  the  impact  of  linearity  on  othe'’ 


types  of  physical  scale  models  deserves  some  clarification. 


The  Maxwell  equations  are  linear  partial  differential  equations.  Since  they 
are  linear  in  the  fields,  these  only  appear  to  the  first  power,  and  linear 
superposition  holds.  Therefore,  any  sum  of  solutions  of  Maxwell's  equations  is 
also  a  solution. 

However,  nonlinear  effects  may  intrude  in  physical  scale  modeling,  and  it  is 
helpful  to  distinguish  between  the  trivial  linearity  of  equations,  the  more 
complicated  linear  response  of  electrical  devices,  and  the  still  subtler  linearity 
of  materials.  Common  examples  of  transitions  from  the  linear  to  the  nonlinear 
regime  occur  in  transducers  used  to  couple  telephone  conversations  with  a  microwave 
beam,  in  crystals  responding  to  intense  laser  beams,  and  in  magnetic  materials. 

Tne  constitutive  relations.  Equations  (5),  permit  a  simple  definition  of  the 
linear  response  of  a  medium.  Accordingly,  materials  are  termed  linear  when  an 
applied  E  or  H  field  induces  a  P  or  M  polarization,  respectively,  that  is 
proportional  to  the  magnitude  of  the  applied  field.  Thus,  our  treatment  also  does 
not  apply  to  substances  where  P  or  M  is  not  proportional  to  the  magnitude  of  the 
applied  field,  such  as  ferroelectrics  or  ferromagnets ,  which  have  nonzero  P  or  M  in 
the  absence  of  any  applied  fields. 

The  Maxwell  equations  are  valid  for  nonhomogeneous  as  well  as  for  homogeneous 
materials.  For'  weak  enough  fields,  however,  different  Carttsian  components  of  E  and 
D  snd  of  B  and  H  are  coupl'''d  through  the  ‘-^lectric  penmi  tt  i  vi  ty  or  dielectric 
tensor,  and  through  the  magnetic  permeability  tensor,  respectively.  These  tensors 
summarize  the  linear  response  of  the  medium  and  are  dependent  on  molecular  and, 
perhaps,  crystalline  structure  of  the  mat^^rial,  as  well  as  on  such  bulk  properties 
as  density  and  temperatur’'’ . 

For  simple  materials,  this  linear  response  may  be  isotropic  in  space, 
whereupon  the  tensors  become  diagonal  with  all  three  elements  equal.  Equations  (5) 


then  become 


^  “  EflE  +  Xg£oE  =  Eod  =  iCg^oE  =  eE 

B  =  yo(H  +  XtnH)  =  Pod  Xn,)H  =  iCmPoH  =  pH 

In  addition,  assuming  no  Hall  or  low  temperature  effects,  Equation  (4)  becomes 
Ohm's  law 


J  =  oE 


where 


Xg  =  electric  susceptibility  (dimensionless) 

K-'e  =  electric  permittivity  relative  to  free  space  (dimensionless) 
e  =  permittivity  of  medium  (F/m) 

Xn,  =  magnetic  susceptibility  (dimensionless) 

<n,  =  magnetic  permeability  relative  to  free  space  (dimensionless) 
(j  =  permeability  of  medium  (H/m) 

0  =  electric  conductivity  of  medium  (mho/m) 


NONLOCALITY  IN  SPACE  AND  TIME 

In  general,  the  basic  linear  connection  between  D  and  E,  as  well  as  between 
H  and  B  and  between  J  and  E,  may  be  nonlocal.  Nonlocality  in  time  and  space 
implies  that  a  principal  property  of  the  relations  (6)  and  (7),  namely  the 
dependence  of  J,  D,  and  B  only  on  values  of  E  and  H  at  the  instant  or  place 
considered,  is  no  longer  true. 

In  general,  the  values  of  J,  D,  and  3  at  a  given  instant  are  determined  not 
only  by  the  values  of  E  and  H  at  that  instant,  but  also  by  values  of  E  and  H  at 
every  previous  instant.  Furthermore,  with  sufficiently  rapidly  varying  applied 
fields,  because  of  the  inertia  of  electrons  and  other  more  massive  atomic  and 
molecular  solid  state  constituents,  establishment  of  P  or  M  polarization  in 
materials  and  the  flow  of  currents  through  them  will  eventually  not  be  abl’-^  to  keep 
up  with  the  changes  in  E  and  H  causing  these  effects  [19-21]. 

Also,  long-range  effects  may  become  important,  so  a  spatially  local  form  of 
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Ohm's  law  is  no  longer  adequate.  Thus,  Equations  (6)  and  (7)  should  be  understood 
as  holding  for  Fourier  transforms  of  the  field  components  in  space  and  time.  The 
permittivity,  permeability,  and  conductivity  tensors  are  then  functions  of 
frequency  f  and  wavelength  A  [19,20].  Nonlocality  in  space  may  often  be  neglected 
in  conductors  up  through  optical  frequencies,  as  long  as  the  mean  free  path  for 
conduction  electrons  colliding  with  the  lattice  structure  is  small  compared  to  the 
skin  depth,  whereupon  e,  Ui  and  a  become  functions  only  of  f  in  dispersive  targets. 
However,  this  f  dependence  can  be  present  down  through  radar  frequencies  and  it 
complicates  scaling,  since  e,  u,  and  o  are  complex  variables  having  both  real  and 
imaginary  (i.e.,  reflective  and  absorptive)  parts. 

ENTRANCE  OF  THE  SKIN  DEPTH 

When  an  electromagnetic  wave  enters  a  good  conductor,  its  amplitude  is  a 
rapidly  decreasing  function  of  distance.  It  is  damped  to  0.36788  (the  reciprocal 
of  the  2.71828  base  for  natural  logarithms)  of  its  initial  amplitude  in  the  skin 
depth  6,  where:  [13i  19] 

1 

6  = -  (8a) 

/  TTf  op 

0  is  also  a  function  of  f,  so  6  decreases  as  either  f  or  o  increases, 
approaching  zero  either  when  o  becomes  very  large  (perfect  conductor)  or  when  f 
becomes  very  large.  A  conductor  with  vanishing  <5  would  reflect  all  incident  light, 
and  would  not  permit  penetration  of  an  electromagnetic  wave  to  any  depth  at,  all. 

We  can  understand  the  electromagnetic  fields  in  the  neighborhood  of  this  thin 
surface  layer,  as  well  as  the  power  absorbed  in  and  scattered  from  it,  by 
considering  waves  incident  on  both  perfect  and  good,  but  not  perfect,  conductors. 

With  waves  incident  on  a  perfect  conductor,  the  boundary  condition  associated 
with  Equation  (2b)  requires  that  surface  charges  move  instantly  in  response  to 
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normal  D  outside  to  produce  the  correct  surface-charge  density  which  gives  zero  E 
inside.  Similarly,  in  response  to  tangential  H  outside,  the  second  boundary 
condition  associated  with  Equation  (2a)  requires  that  surface  charges  move  fast 
enough  to  produce  the  correct  surface  current  to  give  zero  H  inside. 

The  boundary  condition  on  normal  B  associated  with  Equation  (1b),  as  well  as 
the  other  boundary  condition  on  tangential  E  associated  with  Equation  (la),  imply 
continuity  of  these  two  across  the  surface.  Thus,  their  vanishing  values  inside 
require  that  these  two  both  vanish  outside  the  surface.  Therefore,  only  normal  E 
and  tangential  H  can  exist  just  outside  a  perfect  conductor,  and  all  fields  drop 
abruptly  to  zero  inside. 

For  a  good,  but  not  perfect,  conductor  the  surface  fields  and  charges  behave 
approximately  the  same  as  for  a  perfect  conductor  [19,  21].  However,  with  finite 
0,  charges  and  currents  do  not  respond  to  fields  instantaneously,  so  there  is  a 
thin  transitional  surface  layer  in  which  tangential  H  outside  undergoes  the  rapid 
exponential  decay  inside  that  is  measured  by  6.  There  is  also  a  normal  E  outside, 
which  undergoes  an  exponential  decay  inside  similar  to  tangential  H,  but  whose 
energy  density  inside  is  much  smaller  than  the  magnetic  energy  density. 

In  addition,  there  are  small  components  of  E  and  H,  tangential  and  normal  to 
the  surface,  respectively,  outside  the  conductor.  These  also  decrease 
exponentially  inside,  and  result  in  a  power  flow  into  the  conductor,  which  can  be 
calculated  from  Poynting's  vector,  S; 

S  =  E  X  H  (8b) 

which  is  the  power  flow  across  a  unit  area  normal  to  S,  or  the  energy  flux  density. 

Vany  such  quantities  can  be  shown  to  be  proportional  to  6.  For  example, 
Panofsky  and  Phillips  have  shown  [22]  that  the  ratios:  (a)  of  tangential  E  at  the 
metal's  surface  to  its  normal  component  there,  as  well  as  (b)  of  normal  H  to  its 
tangential  component  there,  are  both  approximately  6/A. 
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In  addition,  the  time-average  power  absorbed  (P)  per  unit  area  (a)  equals  the 
time-average  normal  component  of  Poynting's  vector  [19],  which  can  be  written 
dP  uuf6  _ 

—  =  -  |hJ^  (8c) 

da  2 

and  thus  is  proportional  to  6.  This  can  also  be  calculated  from  ohmic  losses  due 
to  energy  dissipation  by  the  surface  currents  in  the  conductor's  transition  layer. 
The  power  absorbed  in  the  conductor  leads  to  a  loss  of  energy  in  the  reflected  wave 
and  a  decreased  scattering  cross  section. 

Therefore,  very  little  power  flows  into  the  conductor,  to  be  dissipated  there 
because  of  wave  damping  due  to  development  of  heat  in  the  surface  resistance. 
Poynting's  vector  at  the  surface  is  very  small,  since  the  ratio  of  E  to  H  within 
the  metal  is  much  smaller  than  its  value  in  free  space,  representing  a  small  energy 
flow  into  the  metal.  The  situation  is  almost  like  that  with  a  perfect  conductor, 
for  which  the  E  vectors  of  the  incident  and  reflected  waves  are  equal  in  magnitude, 
but  exactly  cancel  at  the  metal's  surface,  whereas  the  H  vectors  are  equal  in 
magnitude,  because  of  the  perfect  reflection,  but  add. 

For  a  realistic  material  such  as  copper,  the  above  formula  f or  6  becomes 

6.6  66,000 

5  =  -  cm  =  — ' "  -  urn  (8d) 

/  f  /  f 

yielding  S  =  0.66  ym  at  f  =  10  GHz,  where  the  X  =  2  cm,  and  6  =  0.047  um  at 
f  =  2000  GHz,  where  A  =  1 50  urn.  This  rapid  attenuation  implies  that  current  flows 
in  a  rather  thin  surface  layer  at  high  frequencies,  and  that  6  is  a  very  small 
fraction  of  A.  Thus,  no  complicated  constructive  or  destructive  interference 
occurs  in  this  surface  layer,  and  at  the  metal's  surface  both  tangential  E  and 
normal  H  are  small  compared  to  normal  E  and  tangential  H,  respectively,  according 
to  the  paragraph  above  Equation  (Be). 


The  above  equation  will  also  yield  the  frequency  at  which  6  and  X  eventually 
become  equal: 

c  10“ 

f  =  -  -  -  GHz  =  20.66  GGHz  (8e) 

X  14. an 

where  c  is  the  velocity  of  light  {3*10®  m/sec),  and  1  GGHz  =  10®  GHz.  This 
frequency  corresponds  to 

6  =  X  =  0.00001452  ytn  =  0.1452  X  (8f) 

or  about  one-quarter  the  Bohr  radius  of  the  hydrogen  atom.  Both  6  and  X  appear  to 
be  in  the  quantum  mechanical  regime,  outside  the  range  of  classical  electromagnetic 
phenomena . 

Although  the  above  conclusions  are  based  on  Equation  (8a),  derived  under  the 
assumption  eio  <<  o,  where  u  is  the  angular  frequency  (=  2:if),  dependence  on  X 
becomes  important  in  physical  scale  models  when  the  mean  free  path  of  electrons 
colliding  with  the  conductor's  solid  structure  becomes  comparable  to  5,  which 
usually  does  not  occur  until  one  reaches  optical  frequencies.  Then,  a  spatially 
local  form  of  Ohm's  law  is  no  longer  adequate,  and  long-range  effects  occur  which 
become  important  in  understanding  the  surface  impedance  of  metals. 

A  metal's  surface  impedance  is  the  ratio  of  tangential  E  to  tangential  H  at 
its  surface.  Although  its  value  can  be  calculated  from  classical  electromagnetic 
theory,  at  sufficiently  low  temperatures  the  spatial  nonuniformity  of  the 
electromagnetic  fields  becomes  important.  Consequently,  a  macroscopic  description 
of  a  metal  in  terms  of  e,  y,  and  o  is  no  longer  possible. 

More  precisely,  as  a  metal's  temperature  falls,  its  o  increases,  and  an 
el'ctromagnet ic  field's  depth  of  penetration  into  a  metal  decreases  in  accordance 
with  the  above  expressions  for  6.  At  the  same  time,  the  amplitude  of  lattice 
vibrations  decr’ases  with  a  decrease  in  temperature,  so  the  electronic  collisional 
mean  free  path  increases  and  may  exceed  6  by  a  large  factor.  When  the  mean  free 


path  is  comparable  with  or  much  greater  than  6,  electrons  will  traverse  regions 
with  different  values  of  E  between  collisions,  and  their  resultant  velocity  will 
depend  on  the  changing  values  of  E  along  the  entire  path.  Since  E  varies  over  the 
electron's  path,  o  is  not  constant  over  all  parts  of  the  metal,  so  Ohm's  law  must 
be  replaced  by  a  more  detailed  expression  where  the  spatial  variation  of  E  is  taken 
into  account. 

The  above  behavior  implies  not  only  that  the  current  at  a  point  depends  on  the 
E  field  at  other  points,  but  also  that  the  effects  of  a  boundary  can  be  felt  at 
substantial  depths.  This  leads  to  a  nonexponential  decay  of  the  field  as  it 
penetrates  into  the  medium.  These  and  other  departures  from  the  behavior  described 
by  the  formulas  given  above  for  6  are  known  collectively  as  the  anomalous  skin 
effect  [19''2l],  and  can  be  utilized  to  map  out  the  energy  levels  allowed  to 
electrons  in  metals. 

In  metals,  electrons  that  have  th'^  most  available  energy,  lie  on  what  are 
called  Fermi  surfaces.  For  most  electronic  properties  of  metals,  the  behavior  is 
determined  by  electron  states  very  near  the  Fermi  surface.  This  is  particularly 
true  with  transport  properties,  for  only  electrons  near  the  Fermi  energy  can  find 
an  unoccupied  state  at  a  nearby  energy.  Correspondingly,  with  small  applied  fields 
it  is  only  states  near  this  surface  for  which  the  occupation  number  changes.  Thus, 
energy  bands  in  the  neigiiborhood  of  the  Fermi  surface  are  the  ones  of  primary 
importance. 

To  summarize,  as  f  increases  the  classical  skin  depth  decreases.  When  6 
becomes  comparable  to  the  mean  free  path  of  carriers,  we  reach  the  nonlocal  regime 
described  above,  which  usually  can  be  ignored  in  conductors  up  through  optical 
frequencies.  When  6  is  much  smaller  than  the  carrier  mean  free  path,  we  encounter 
what  is  called  the  extreme  anomalous  skin  effect.  More  detailed  study  of  the 
anomalous  skin  effect  leads  to  information  about  the  transition  between  the 


classical  skin  effect  and  the  extreme  anomalous  region,  and  to  information  about 
details  of  the  shape  of  the  Fermi  surface  [20]. 

NONLINEAR  OPTICS 

Conductors  are  composed  of  atoms  and  molecules  bound  together  in  a  regular 
arrangement  to  form  a  solid.  The  electrons  in  the  outer  shells  of  these  atoms  and 
molecules  can  be  regarded  as  being  in  a  potential  in  which  both  the  electrons  and 
the  resulting  ions  will  oscillate  linearly  at  small  amplitudes. 

Materials  showing  a  linear  response  to  weak  electromagnetic  fields  will 
eventually  show  nonlinear  behavior  at  higher  field  strengths  when  electronic  or 
ionic  oscillators  are  driven  to  large  amplitudes.  A  monochromatic  wave,  when 
pulsed,  develops  a  spectrum  of  frequencies  which,  for  large  amplitude,  may  then 
generate  waves  in  the  medium  with  harmonic,  sum  and  difference,  cubic,  and  higher 
nonlinear  terms,  giving  rise  to  a  broader  spectrum  of  frequencies.  With  the 
development  of  lasers,  such  nonlinear  behavior  has  become  an  active  research  area, 
called  nonlinear  optics  [19,  23]. 

Lasers  are  capable  of  generating  pulses  with  peak  E  fields  of  10  or  even  100 
Gigavolts/cm.  On  the  other  hand,  the  static  E  field  felt  by  an  electron  in  the 
hydrogen  atom's  ground  state  is  only  5.1^  GV/cm,  and  fields  of  100  kV/cm  are  equal 
to  the  average  internal  fields  seen  inside  a  dielectric  lattice.  Since  fields  of 
100  kV/cm  can  thus  cause  electrical  breakdown  in  a  solid  and  destroy  it,  the  laser 
fields  often  used  with  physical  scale  models  may  well  be  capable  of  driving  atomic 
oscillators  into  their  nonlinear  range. 

MODELING  OF  ELECTROMAGNETIC  SYSTEMS 

Scaling  laws  result  from  imposing  exact  conditions  on  the  Maxwell  equations. 

If  a  ni-’d'^l  is  to  be  an  accurate  simulation  of  a  full-scale  system,  transformations 
relating  model  to  full-scale  quantities  will  transform  Maxwell's  equations  from  one 
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system  to  another.  Since  electromagnetic  fields  in  both  systems  must  satisfy 
Maxwell's  equations,  such  a  transformation  will  determine  how  e,  u,  and  o  must 
transform  for  simulation  to  be  accurate.  Replication  of  each  coefficient  of  the 
Maxwell  equations,  considered  separately,  will  then  give  three  general  nonlinear 
modeling  equations  in  six  variables. 

To  obtain  these  equations,  consider  a  linear,  homogeneous,  isotropic  system 
with  no  free  charge,  where  Equations  (6)  and  (7)  are  valid,  and  Maxwell's  two  curl 
equations,  given  in  Equations  (1a)  and  (2a),  take  the  form 

3H 

V  X  E  =  -  u  —  {9a) 

at 


aE 

VxH=oE+e—  (9b) 

at 


The  electromagnetic  properties  of  a  configuration  of  imperfect  conductors, 
dielectrics,  and  magnetic  materials  in  unbounded  space  can  be  compared  with  the 
properties  of  a  scale  model,  which  differs  geometrically  from  the  original  only  in 
having  its  dimensions  changed  by  a  factor  i,  which  may  be  larger  or  smaller  than 
one.  If  model  characteristics  are  denoted  by  primes,  each  model  length  L'  is 
related  to  a  corresponding  length  L  in  the  original  system  according  to 

L'  =  IL  (10a) 

The  conditions  imposed  by  Equation  (10a)  represent  the  requirements  for  a 
mechanical  model  of  the  material  portions  of  a  full-scale  system,  in  which  there  is 
geometrical  similarity  in  shapes  of  corresponding  material  parts.  Frequently,  a 
quantity  p  is  defined,  the  ratio  of  any  full-scale  to  corresponding  model  length, 
called  the  mechanical  scale  factor.  The  parameter  p  is  normally  chosen  to  yield  a 
model  of  convenient  size,  would  usually  be  greater  than  one  for  smaller  models  used 
at  higher  f,  and  satisfies 

p  =  1  /«,  (10b) 
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Time  is  likewise  scaled 


1 

t'=Tt=-t  (10c) 

f 

where  T  is  the  ratio  of  model  characteristic  times  to  those  in  original  system  and, 
since  the  field's  period  is  inversely  proportional  to  frequency,  we  redefine  f  to 
be  the  ratio  of  model  to  original  frequencies. 

The  requirement  that  a  model's  electromagnetic  properties  shall  be  the  same  as 
those  of  the  original,  except  for  a  change  of  scale,  means  that  the  E  and  H  fields 
in  the  space  about  a  model  will  pick  up  multiplicative  factors 

E'  =  e  E  and  H'  =  h  H  { 1 1  a) 

where  e  and  h  are  the  ratio  of  model  characteristic  electric  and  magnetic  field 
intensities  to  those  in  original  medium,  respectively.  Likewise,  c,  m,  and  o  have 
new  values  in  the  model  denoted  by  primes. 

Since  Equation  (10c)  states  that  all  times  are  to  be  multiplied  by  a  factor  T 
in  a  model,  first  time  derivatives  may  be  replaced  by 

3  13  3 

—  - - -  f  —  (11b) 

3t'  T  3t  3t 

The  curl  operator,  which  involves  first  derivatives  with  respect  to  space 
coordinates,  may  similarly  be  replaced  according  to  Equation  (10a).  Then,  starting 
with  model  variables,  the  Maxwell  Equations  (9)  take  the  form 

_  _  1  _  _  e  _  _  3H' 

V  X  E'  =  -  V  X  E'  =  -  V  X  E  =  -  u’  — 

I  I  3t' 

(12a) 

1  _  _  n  _  _  _  3E' 

V  X  H'  =  -  7  X  H'  =  -  V  X  H  =  o'S'  +  e'  — 

I  I  3t' 

Multiplying  the  first  equation  by  (H/e),  the  second  by  (il/h),  and  continuing 
the  transformation  from  primed  to  unprimed  variables 


h  3H 

-  u'if  - 

e  3t 


(12b) 


I  3H' 

VxE=--u'  —  = 
e  3t' 

_  _  il  _  3E'  e  _  e  3E 

V  X  H  =  -{o'E'  +  e'  —  }  =  oH  -  E  +  e’if - 

h  3t'  h  h  3t 

Since  the  electromagnetic  boundary-value  problems  presented  by  the  model  and 
the  original  are  to  be  similar,  we  can  determine  how  e,  u,  and  o  must  transform  by 
requiring  that  the  coefficients  of 

3H  _  3E 

—  ,  E  ,  and  — 

3t  3t 

be  identical  in  both  model  and  original. 

Equating  the  coefficients  in  Equations  (12b)  to  those  in  Equations  (9),  three 
obviously  independent  relations  result 

h 

M  =  u’if  - 
e 

e 

0=0'?,  -  (13a) 

h 

e 

E  =  E ' f  - 

h 

For  simplicity  of  notation,  since  I,  f,  e,  and  h  are  ratios  of  model  to 
original  sines,  frequencies,  E  and  :i  fi-^lls,  respectively,  we  now  redefine  e,  u, 
and  a  t,o  be  ratios  of  model  permittivity,  permeability,  and  conductivity  to  those 
of  the  original. 

I’hen,  the  above  thret?  equal  iem'  t  t!ie  simpler  form 


and  involve  the  seven  variables 


e 

c ,  u ,  a,  I,  f ,  and  - 
h 

Since  the  E  and  H  field  coefficients  only  appear  as  a  ratio,  we  then  define 
0 

C  =  -  (U) 

n 

which  is  actually  the  ratio  of  model  impedances  to  those  in  the  original,  and  our 
thr-''e  '■’quit  ions  take  the  final  form 

M’-f  =  C  (15a) 

ok  =  1  (15b) 

cUc  =  1  (15c) 

in  a  total  of  six  variables. 

SOi.UTIONS  OF  THE  MODEMNG  EQUATIONS 

Complete  solutions  to  the  three  nonlinear  equations  in  six  v  iriables  given  in 
Equa''ion;;  i’15)  have  not  apparently  been  published.  In  practice,  one  seeks 
■se’etions  for  I,  f,  and  r,  as  functions  of  the  parameters  e ,  u,  and  o  available  in 
ar'tual  materials  used  for  models. 

nK'lKHAl.  SOLUTION 

Eliminating  t,  the  ratio  of  fields,  from  Equations  (15)  yields  the  three 
express i ons 


The  first  two  are  similar  to  Stratton's  conditions  for  electrodynamic 
similitude  [133.  although  arbitrary  constants  have  been  evaluated  and  no  appeal  is 
made  to  dimensionless  measure  numbers  of  the  field  variables. 

The  third  equation  is  new,  although  it  has  often  been  stated  that  scaled  model 
conductivity  should  be  p  times  that  in  the  original  system  [1,  4,  5,  8-11,  14-16, 


What  has  been  untouched  is  other  solutions  of  Equations  (15)  and  the  behavior 
of  electromagnetic  fields  in  physical  scale  models.  Going  back  to  Equations  (15), 


and  solving  each  of  these  expressions  for  ^ 


(17a) 


(17b) 


(17c) 


However,  the  presence  of  all  six  variables  in  Equations  (17)  complicates 
interpretation. 

On  the  other  hand,  if  is  eliminated  from  Equations  (15),  two  new  expressions 


result 


(I8a) 


(I8b) 


A  third  is  also  produced  which  is  actually  the  same  as  Equation  (16c). 


►  AA 


If  f,  last  of  the  three  variables,  is  eliminated  from  Equations  (15),  only  one 
relation  is  produced,  which  is  the  same  as  Equation  (18a). 

This  constitutes  a  complete  set  of  solutions  to  the  three  nonlinear  Equations 
(15),  as  can  be  seen  by  successively  eliminating  either  I,  or  f,  or  ^  from 
Equations  (15)  and  finding  all  solutions  for  each  of  the  remaining  two  variables  as 
a  function  of  the  third  and  of  the  remaining  three  parameters  e,  u,  and  o. 

Equations  (18)  may  have  profound  implications  for  physical  scale  models, 
because  they  provide  a  general  solution  relating  the  ratio  of  electric  and 
magnetic  fields  as  defined  in  Equation  (1*1),  to  f  and  the  three  material  parameters 
e,  w,  and  a.  Equations  (16),  (17),  and  (18)  will  now  be  used  to  show  how  finite  a 
affects  the  surface  currents  and  to  discuss  proper  scaling  of  the  electromagnetic 
fields  and  the  scattering  cross  section,  differences  between  geometric  and  complete 
scaling,  and  dependence  of  absorption  and  energy  balance  on  f. 


FINITE  CONDUCTIVITY  AND  SURFACE  CURRENTS 

In  practice,  scale  models  are  limited  by  values  of  e,  u,  and  o  obtainable  in 
model  materials.  For  many  measurements,  e'  and  y'  in  the  model  will  not  differ 
appreciably  from  their  values  in  the  full-scale  system.  According  to  Equation 
(I8a),  the  ratio  of  fields  should  also  not  differ  appreciably  in  a  completely 
scaled  system. 

However,  if  one  uses  Equation  (1*)),  definition  of  the  impedance,  in  Equation 
(I8b) ,  one  obtains 


=  e 


(19) 


0 

^  V  uf 

which  implies  that,  if  o  (or,  more  properly,  o/y)  does  not  increase  linearly  with 
f,  the  ratio  of  H  fields  shrinks  with  respect  to  the  ratio  of  E  fields. 


This  is  understandable  because  a  a  not  increasing  fast  enough  for  proper 


scaling  leads  to  diminished  surface  currents,  which  cause  a  decline  in  their 


associated  H  fields. 

The  coupling  between  electromagnetic  fields  and  surface  currents  is  simple  and 
direct . 

For  a  perfect  conductor,  tangential  E  at  the  surface  is  zero  (otherwise,  any 
surface  current  would  be  infinite),  but  tangential  H  is  finite,  and  the  fields  are 
damped  waves  that  fall  off  infinitely  rapidly  with  penetration  into  the  metal,  so  H 
and  E  are  both  zero  directly  below  the  surface.  Stokes's  theorem  implies  a 
surface-current  density,  numerically  equal  to  tangential  H,  but  at  right  angles  to 
it,  to  account  for  the  rapid  decrease  of  tangential  H  from  its  finite  surface  value 
to  zero  directly  below  the  surface.  This  current  can  flow  without  a  corresponding 
tangential  E  because  of  the  perfect  conductivity. 

Inside  a  good  conductor,  the  transverse  H  falls  off  rapidly  as  we  penetrate 
the  metal.  This  rapid  change  of  H  with  depth  leads  to  a  large  V  x  H,  and  hence  the 
Maxwell  Equation  (2a)  implies  a  large  current  density  J,  since  conduction  currents 
are  much  larger  than  displacement  currents  in  a  good  conductor.  Thus,  there  will 
be  a  large  current  flowing  in  a  thin  metal  surface  layer.  By  Ohm's  law  [Equation 
(7)],  this  surface  current  is  parallel  to  E  and  proportional  to  it,  but  the 
magnitude  of  E  is  small  for  large  values  of  o.  Since  H  drops  to  zero  well  inside 
the  conductor,  Stokes's  law  implies  the  tangential  component  of  H  just  outside  the 
metal  surface  equals  the  current  flowing  parallel  to  the  surface,  per  unit  length 
of  surface. 

In  the  limit  of  infinite  conductivity,  J  in  the  surface  layer  becomes 
infinite,  but  5  becomes  zero  in  such  a  way  that  a  finite  current  per  unit  length  of 
surface  results,  which  is  numerically  equal  to  the  tangential  component  of  H  at  the 
surface.  On  the  other  hand,  tangential  E  goes  to  zero  all  through  the  surface 
layer,  so  tangential  E  at  the  surface  is  also  zero. 
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In  fact,  when  both  target  and  model  conductivity  are  infinite,  so  both  metal 
surfaces  are  perfect  reflectors  for  waves  at  all  frequencies,  then  a  becomes  the 
ratio  of  two  infinities,  which  is  Indeterminate.  Therefore,  one  general  scaling 
equation.  Equation  (15b),  becomes  trivially  satisfied,  and  the  other  two.  Equations 
(15a)  and  (15c),  have  the  simple  solutions  given  in  Equations  (l6a),  (17a),  (17c), 
and  (l8a).  Furthermore,  since  I  and  f  appear  only  in  the  product  (Hf)  in  the 
general  modeling  Equations  (15a)  and  (15c),  if  one  assumes  geometric  scaling 

if  =  1  (20) 

then,  according  to  Equation  (l8a),  the  usual  f-independence  assumed  for  e  and  y 
will  imply  electromagnetic  fields  whose  ratio  c  doesn't  depend  on  f. 

In  addition,  with  Equation  (20)  used  to  link  I  and  f  in  the  general  modeling 
Equations  (15a)  and  (15c),  all  dependence  on  both  these  variables  disappears  from 
their  solutions,  so  no  f-dependence  of  material  properties  is  required  for  proper 
scaling. 

On  the  other  hand,  with  both  target  and  model  a  taken  to  be  finite,  additional 
requirements  imposed  by  the  general  scaling  equation.  Equation  (15b),  lead  to  a 
patent  f-dependence  of  material  properties  being  needed  for  proper  scaling  that 
complicates  the  new  solutions  given  in  Equations  (I5b),  (I6c),  (17b),  (I8b),  and 

(19) .  To  better  understand  this  required  f-dependence  of  material  properties,  it 
is  helpful  to  consider  the  differences  between  geometric  and  complete  scaling. 

DIFFERENCES  BETWEEN  GEOMETRIC  AND  COMPLETE  SCALING 

In  our  notation,  the  geometric  scaling  condition  is  expressed  by  Equation 

(20) .  Increased  f  is  then  associated  with  decreased  size,  so  a  scaled  target  has 

the  same  dimensions  in  scaled  wavelengths  as  the  actual  target  has  in  actual 
wavelengths . 

For  realistic  nat^rials  with  finite  o,  models  built  according  to  the  geometric 


scaling  condition  of  Equation  (20)  lead  to  further  conditions  on  ini'- 
properties  involved  if  complete  scaling  is  to  prevail.  Thf'Si'  von  d  i ’■  i  ons  ■ -.r,  ; 
reatly  obtained  by  using  Equation  (20)  to  link  '<  and  f  in  the  gener-al  mcxv-.inr 
Equations  (15).  This  process  permits  elimination  of  1.  Solutions  for  ;  :n  tv 
be  obtained  as  functions  of  f  and  of  the  parameters  c,  u,  and  j  available  in  aetijal 
model  materials. 

Explicitly,  using  Equation  (20)  to  eliminate  I  from  Equations  (15),  the  three 
independent  basic  modeling  relations  simplify  to 
5  =  u 

=  f  (21) 

=  1 

Eliminating  c,  from  Equations  (21)  yields  the  analogues  of  Equations  (16) 

UE  =  1 

ou  =  f  (22) 

0 

cf 

These  are  the  conditions  Maxwell's  equations  place  on  the  materials  needed  for 
accurate  simulation  of  systems  by  models.  We  see  that  various  combinations  of  o 
with  e  and  p  should  Increase  linearly  with  f,  while  c  and  u  also  satisfy  the  first 
of  Equations  ( 22) . 

In  Equations  (21),  no  I  or  f,  common  to  more  than  one  equation,  can  be 
eliminated  there.  However,  solving  each  of  Equations  (21)  for  c,  the  analogues  of 
Equations  (17)  are 


and  multiplying  pairs  of  these  yields  the  analogues  of  Equations  (18) 


Equations  (24)  are  very  similar  to  Equations  (18),  and  again  imply  that  a  o 
not  increasing  fast  enough  for  proper  scaling  leads  to  diminished  surface  currents 
and  declining  H  fields. 


CONDITIONS  FOR  AN  ABSOLUTE  MODEL 

Sinclair  [5]  also  discussed  modeling  theory  for  electromagnet ic  systems. 

Since  there  are  only  four  fundamental  units  (mass,  length,  time,  and  charge),  four 
scale  factors  should  suffice  to  describe  any  electromagnetic  quantity.  Thus,  he 
introduced  four  factors  to  relate  full-scale  length,  frequency  (or  time),  E,  and  H, 
to  model  quantities,  analogous  to  those  in  Equations  (10)  and  (11). 

He  derived  three  equations  similar  to  Equations  (13a),  intending  to  absolutely 
replicate  a  full-scale  system,  with  both  geometries  and  field  strengths  being 
modeled.  He  then  calculated  how  24  quantities  [including  Poynting  vector  and  radar 
cross  section  (RCS)]  transform  as  functions  of  the  four  scale  factors,  and  claimed 


h*^  could  obtain  quantitative  data  on  all  electromagnetic  properties  of  a  system, 


including  absolute  power  levels  as  well  as  configurations  of  lines  of  force. 

However,  Maxwell's  equations  fix  only  the  ratio  c  =  e/h,  and  not  e  or  h, 
separately.  Thus,  Sinclair  simply  noted  that  if  a  specific  value  is  assigned  to 
either  e  or  h,  his  model  is  absolute  and  will  yield  quantitative  results  for  all 
quantities.  Otherwise,  his  model  is  a  geometric  one. 

Many  scale  models  are  built  to  satisfy  the  geometric  scaling  condition  of 
Equation  (20),  so  one  cannot  then  deduce  power  levels  directly  from  model 
measurements.  However,  absolute  values  can  be  obtained  by  measuring  relative  power 
in  the  model  system  and  calibrating  this  power  by  comparing  signals  returned  from  a 
model  with  those  from  a  target  of  known  RCS. 

Unfortunately,  Sinclair  [5l  provided  no  general  solutions  of  the  modeling 
equations,  and  did  not  consider  effects  of  differences  between  geometric  and 
complete  scaling  on  electromagnetic  fields,  nor  the  influence  of  finite  o  on 
surface  currents.  Also,  little  attention  was  paid  to  dependence  on  frequency. 

ABSORPTION  IN  NONDISPERSIVE  TARGETS 

The  scaling  laws  considered  above  imply  replication  of  each  field  coefficient 
in  the  Maxwell  Equations  (9),  considered  separately,  and  consist  of  three  general 
nonlinear  modeling  Equations  (15)  in  the  six  ratios:  I,  f,  c,  e,  u,  and  a. 

However,  now  going  back  to  the  original  meaning  of  e,  u,  and  a  (and  f),  while  their 
static  values  are  reasonably  well  known,  absorption  at  higher  frequencies  occurs 
near  eigenfrequencies  of  the  molecular  or  electronic  vibrations  causing  P  or  M 
polarization  of  materials. 

At  these  frequencies,  the  linear  constitutive  Equations  (6)  and  Ohm's  law, 
Equation  (7),  are  not  valid.  The  values  of  these  eigenfrequencies  depend  on  the 
substance  concerned,  and  vary  widely,  but  become  important  whenever  e,  u,  and  o 
exhibit  the  f-dependence  known  as  dispersion.  Moreover,  since  the  presence  of 


dispersion  in  general  signifies  a  dissipation  of  energy,  nearness  to  these 
eigenfrequencies  is  also  indicated  whenever  materials  absorb.  Thus,  a  dispersive 
medium  is  also  an  absorbing  medium.  In  addition,  the  frequencies  at  which 
dispersion  phenomena  first  appear  may  be  entirely  different  for  the  electric  and 
magnetic  properties  of  a  substance,  as  well  as  being  dependent  on  direction  in  the 
lattice.  In  this  report,  while  we  continue  to  acknowledge  anisotropy,  we  defer  its 
proper  treatment  until  a  future  time. 

Since  absorption  at  higher  frequencies  occurs  near  eigenfrequencies  of  the 
molecular  or  electronic  vibrations  causing  P  or  M  polarization  of  materials,  one 
would  expect  this  absorption  to  be  associated  with  \  for  the  electromagnetic  fields 
eventually  becoming  comparable  with  atomic  dimensions,  so  the  macroscopic 
formulation  utilized  above  is  invalid.  On  the  contrary,  there  is  an  extensive 
frequency  range  in  which  absorption  and  dispersion  phenomena  are  important  but, 
fortunately,  X  for  electromagnetic  fields  is  much  larger  than  that  of  atomic  or 
molecular  dimensions,  so  a  macroscopic  formulation  utilizing  e,  p,  and  o  is 
appropriate. 

The  quantum-mechanical  band  structure  of  solids  shows  that  conductors  have 
"free"  electrons  in  partially  filled  bands,  while  Insulators  have  bands  filled  to 
the  full  extent  permitted  by  the  Pauli  principle.  In  conductors,  these  "free" 
electrons  have  their  motion  damped  by  collisions  involving  momentum  transfers  to 
the  atoms  making  up  the  resistive  medium,  lattice  imperfections,  and  impurities. 
This  electron  inertia  usually  is  felt  at  a  lower  frequency  than  any  resonant 
absorption  frequency,  and  is  then  the  first  cause  of  dispersion.  The  Drude  theory 
explains  these  dispersive  effects  of  electron  inertia  by  means  of  a  complex  e  [19, 
29].  Hence,  it  is  quite  proper  to  consider  nondispersi ve  absorption  up  through 
microwave  frequencies  to  the  range  where  the  effects  of  electron  inertia  begin  to 


ELECTRIC  AND  MAGNETIC  SCATTERING 


Since  both  the  physical  causes  and  mathematical  formulation  of  electric  and 
magnetic  scattering  are  quite  different,  each  type  of  scattering  can  be  considered 
separately. 

EJectric  Scattering 

For  nondispersi ve  materials  with  real  e  and  o,  typical  electric  absorption 
would  take  place  in  a  lossy  dielectric,  whose  parameters  satisfy  the  Maxwell  curl 
equation 

3E 

V  X  H  =  oE  +  e  —  (9b) 

3t 

The  two  terms  on  the  right  hand  side  of  Equation  (9b)  represent  conduction  and 
convection  currents,  which  depend  on  frequency.  However,  any  variable  field  can  be 
Fourier  transformed  into  a  series  of  single-frequency  components,  in  which  all 
quantities  depend  on  time  through  the  complex  exponential  factor  exp(iu)t).  This 
form  makes  it  possible  to  reduce  operations  of  time  differentiation  to  simple 
multiplication  by  iw,  and  also  to  more  readily  consider  all  parameters  of  equations 
as  complex  variables. 

Assuming  an  exp(iut)  time  dependence,  monochromatic  fields  are  singled  out, 
and  Equation  (9b)  becomes 

V  X  H  =  oE  +  iewE  (25a) 

whereupon  it  can  be  seen  that  the  conduction  and  convection  currents  are  93  degrees 
out  of  phase.  The  lo  dependence  apparent  on  the  right  hand  side  implies  that  the 
resistive  conduction  current  dominates  at  lower  ui,  whereas  the  reactive  convection 
ciirrent  dominates  at  higher  u. 

If  0  is  considered  to  be  the  seat  (dT  all  electnic  losses,  it  is  convenient  to 
define  a  quantity  n  called  the  ■'electric)  loss  tangent 
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0 

conduction  current  density 

(26) 

WE 

displacement  current  density 

Since  n  is  obviously 

a  function  of  frequency,  then  the  degree  to  which 

a 

material  exhibits  conductor  or  dielectric  properties  depends  on  frequency. 

More 

specifically. 

n  f"' 

conductor 

\  «  ^ 

dielectric 

On  the  other  hand,  scaling  of  o  and  e  according  to  Equations  (13a),  which 
insures  proper  replication  of  each  coefficient  of  the  Maxwell  equations,  leads  to  a 
loss  tangent  not  changing  with  frequency,  so  the  degree  to  which  a  properly  scaled 
material  exhibits  the  properties  of  a  conductor  or  dielectric  should  not  change 
with  frequency. 

Magnetic  Scattering 

Typical  magnetic  absorption  occurs  in  magnetized  ferrites,  which  can  be  very 
anisotropic  (e.g.,  directional  couplers)  and  whose  parameters  may  show  strong 
nonlinearity  in  their  dependence  on  the  magnitude  of  H.  However,  the  right-hand 
side  of  the  other  Maxwell  curl  equation 

3H 

V  X  E  =  -  M  —  (9a) 

at 

consists  of  only  one  term. 

4ssuning  a  time  dependence  of  the  form  expdwt),  monochromatic  fields  are 
again  singled  out,  and  Equation  (9a)  becomes 

V  X  E  =  -  ipwH  (25b) 

whereupon  it  can  be  seen  that  nonriispersive  materials  with  real  u  have  no  room  for 
magnetic  absorption.  However,  magnetic  absorption  will  occur  in  dispersive 


materials  when  u  acquires  an  imaginary  part.  The  imaginary  part  of  u  then  plays  a 


role  In  Equation  (25b)  analogous  to  that  of  o  in  Equation  (25a),  so  it  becomes 


possible  to  define  a  magnetic  loss  tangent  that  is  the  ratio  of  two  current 
densities  like  those  in  Equation  (26). 

However,  for  nondispersive  materials  with  real  p,  there  is  no  magnetic 
conduction  current,  as  well  as  no  free  magnetic  charge,  and  the  magnetic  loss 
tangent  is  zero  at  all  frequencies. 

Furthermore,  scaling  of  u  according  to  Equations  (13a)  does  not  allow  for 
magnetic  absorption  here  at  any  frequency. 

PLANE  WAVES  IN  DISSIPATIVE  MATERIALS 

The  electric  and  magnetic  properties  of  materials  both  influence  the 
reflection  and  absorption  of  electromagnetic  waves  by  targets.  We  now  consider  how 
the  reflection  coefficient  and  radar  cross  section  (RCS)  change  as  frequency 
increases,  and  as  incomplete  scaling  causes  n  to  decrease  through  the  value  n  =  1 
[13.  19,  25]. 

Both  E  and  H  satisfy  (damped)  wave  equations  of  the  same  form.  Taking  the 
curl  of  Equation  (9a),  substituting  from  Equation  (9b)  for  V  x  H,  and  assuming  no 
free  charges  on  surfaces,  yields 

9E  3^E 

V^E  -  op - ep  -  -  0  (27a) 

at  at" 

An  analogous  procedure  yields  an  identical  equation  for  H. 

In  a  rectangular  coordinate  system,  for  incident  uniform  plane  waves  normal  to 
a  surface  at  z  =  0,  the  solution  of  Equation  (27a)  takes  the  form 

E  =  Eo  exp(i(Dt  -  Tz)  (28) 

where  Eq  is  a  constant  amplitude  vector  in  the  direction  of  the  wave's 
polarization,  and  Y  is  the  complex  propagation  constant..  Inserting  Equation  (28) 
into  Equation  f27a)  yields 


=  ioipCo  +  iwe) 


=  -  w*eij(1  -  in)  =  -  k^d  -  in)  (29) 

Thus,  Y  is  actually  a  function  of  only  two  parameters:  (a)  k,  the  propagation 

constant  of  light  when  o  =  0 

_  f  2?: 

k  =  u)/eu  =  2Tr-  =  —  (30) 

V  X 

where  v  is  the  velocity  of  light  in  a  material  medium  (v  =  1//eu),  and  (b)  the 
electric  loss  tangent  n  from  Equation  (26). 

By  taking  the  square  root  of  Equation  (29),  Y  itself  can  be  written  as  a  sum 
of  the  proper  real  (a)  and  imaginary  (8)  parts 

Y  =  a  +  iB  (31 ) 

to  make  Equation  (28)  a  solution  of  Equation  (27a).  The  nonnegative  character  of 
€,  y,  and  a  causes  Y^  to  lie  in  the  second  quadrant  of  the  complex  pxane,  so  its 
positive  square  root  lies  in  the  first  quadrant,  and  thus  a  and  8  have  the  same 
sign.  The  other  square  root  would  be  denoted  by  -Y. 

Substituting  expressions  of  the  form  of  Equation  (28)  for  E  and  H  into 
Maxwell's  equations,  the  plane  wave  considered  is  polarized,  with  E  and  H  at  right 
angles  to  each  other  and  to  the  direction  of  propagation.  The  wave  (or 
characteristic)  impedance  Z  of  a  medium  is  defined  as  the  ratio  of  the  magnitude  of 
E  to  the  magnitude  of  H,  where 
icoy 

Z2  =  - 

0  +  iwe 

y  1  1 

= - =  (Zo)^  -  (32) 

E  (1  -  in)  (1  -  in) 

Thus,  Z  also  turns  out  to  be  a  function  of  two  parameters:  (a)  Zp,  the 
characteristic  impedance  of  a  medium  when  o  =  0 


30 


Zo  = 

which  reduces  to  376.6  fi  for  free  space,  and  (b)  the  electric 
Equation  C26) . 


(33) 

loss  tangent  n  from 


3y  taking  the  square  root  of  Equation  (32),  Z  itself  can  also  be  written  as  a 
sum  of  real  (Z')  and  imaginary  (Z")  parts 

Z  =  Z'  +  iZ”  (3i4) 

where  Z  has  the  dimensions  of  a  resistance.  Since  it,  like  T,  is  complex  in  a 
conducting  medium,  th^re  is  in  general  a  phase  angle  between  the  time  dependence  of 
E  and  H.  The  nonnegative  character  of  e,  u,  and  o  now  causes  Z^  to  lie  in  the 
first  qi'adrant  of  the  complex  plane,  so  its  positive  square  root  lies  in  the  first 
octant,  and  thus  Z'  and  Z"  have  the  same  sign. 

The  properties  of  the  medium  may  be  neatly  expressed  in  terms  of  Y  and  Z, 
allowing  the  two  Fourier-transfot -led  Maxwell  curl  Equations  (25a)  and  (25b)  to  be 
written  _  _  _ 

V  V  (ZH)  =  YE  (35a) 

7  X  E  =  -  Y(ZH)  (35b) 


From  them,  as  wi^'^  Eqtiation  (27a),  we  may  derive  the  wave  equations 

V^E  -  Y^E  =  0  (27b) 

V^H  -  Y^H  =  0  (27c) 

satisfied  by  all  components  of  E  and  il. 

ftll  solutions  of  Equations  (27b),  (27c),  and  (35)  are  thus  functions  of  only  Y 
ind  and  th^’oefore  of  the  three  parameters:  k,  Zg,  and  n,  as  given  by  Equations 
( and  (26).,  respectively.  It  is  often  convenient  to  use  the  Maxwell 


o'lu '.t  i "no,  and  th'^  wave  equations  in  these  forms. 

Sd.Anl'j'l  'T'  THE  .SCATT.ERIM'l  CROSS  SECTION 


iinf'c  all  solutions  of  the  Fourier-transformed  wave  and  Maxwell  equations  for 


E  and  H,  Equations  (27b),  (27c),  and  (35),  are  functions  only  of  Y  and  Z,  and 
therefore  of  the  three  parameters  k,  Zo,  and  n,  investigating  the  scaling  of  any 
quantity  derived  from  the  E  and  H  fields  reduces  to  investigating  the  scaling  of 
these  three  quantities,  given  by  Equations  (30),  (33) i  and  (26),  respectively. 

For  nondispersi ve  materials,  neither  e'  nor  u'  will  differ  from  their  values 
in  the  full-scale  system.  As  we  have  seen  from  Equations  (18)  and  (2A) ,  the  ratio 
of  fields  should  also  not  differ  in  a  completely  scaled  system.  Thus,  Zq  should 
not  change  with  frequency. 

Equations  (22')  are  the  conditions  Maxwell's  equations  place  on  the  materials 
needed  for  accurate  simulation  of  systems  by  models.  If  neither  c’  nor  p'  differ 
from  their  values  in  the  full-scale  system,  the  first  of  Equations  (22)  is 
automatically  satisfied,  and  the  second  and  third  of  Equations  (22)  both  imply  that 
0  s'noiild  increase  linearly  with  f.  Therefore,  neither  n  nor  Z  should  change  with 
frequency  in  a  completely  scaled  system. 

Thus.,  under  complete  scaling,  neither  Z,  Zq,  nor  n  will  change  with  frequency, 
while  Y  and  k  bo*'h  can  easily  be  seen  from  Equations  (29)  and  (30)  to  increase 
linearly  with  f. 

On  the  ether  hand,  for  the  special  case  of  a  uniform  plane  wave  normally 
inci.  dent  from  free  ;:pace  (mediiim  1)  upon  the  plane  surface  of  a  dissipative  block 
(m-'dium  O'),  such  g'^nerality  is  not  necessary. 

In  m'Cdium  l,  E  consists  of  incident  and  reflected  waves,  and  one  defines  a 
refl-'^cticn  co^^fficient  R,  measuring  the  amplitude  ratio  of  reflected  to  incident  E. 
In  m-'di  in  thei^e  is  also  a  transmitted  wave,  so  one  can  define  a  transmission 
coeffioi-nt  T,  analogous  to  R.  Both  Stratton  [13]  and  Adler  et  al .  [25]  calciulated 
R  ini  T  for  scattering  from  dissipative  surfaces  and  found  that  both  these 
co'*f f ici.-^nts  depend  only  on  the  Z's  and  not  on  the  Y's. 

Thus,  under  complete  scaling,  since  both  R  and  T  depend  on  the  Z's  and  not  on 


the  y's,  although  the  scaled  RCS  of  a  target  will  decrease  as  the  square  of 
increasing  frequency  (like  a  physical  area),  its  properly  scaled  reflection 
coefficient  should  not  change. 

On  the  other  hand,  under  incomplete  scaling,  with  the  scatterer's  u,  e,  and  o 
not  regarded  as  functions  of  frequency,  and  with  medium  1  being  free  space  with 
loss  tangent  zero,  the  only  remaining  f-dependence  resides  in  n  for  medium  2,  which 
is  also  a  function  of  the  values  of  a  and  e  there. 

ENERGY  BALANCE  AND  DEPENDENCE  ON  FREQUENCY 

It  is  physically  impossible  to  satisfy  scaling  conditions  on  o  when  good 
conductors  like  copper  or  aluminum  are  used  in  the  original  system  and  the 
frequency  change  is  appreciable  because,  for  copper  and  aluminum,  o  is  so  near  the 
theoretical  maximum  that  it  cannot  be  increased  in  accordance  with  Equations  (I6c) 
and  (I8b).  However,  water,  resistive  conductors,  and  lossy  dielectrics  have 
conductivities  far  enough  from  the  theoretical  maximum  that  their  a  can  be 
increased  in  accordance  with  Equations  (I6c)  and  (l8b)  for  appropriate  changes  in 
f.  Experience  with  more  severe  requirements  in  obtaining  materials  with  necessary 
0  is  found  in  geophysical  applications,  where  p  in  Equation  (10b)  must  be  of  order 
200  to  1000  for  modeling  mineral  prospecting  problems  and  perhaps  as  large  as 
100,000  to  1,000,000  to  study  effects  of  large-scale  structures  on  geomagnetic 
variations  [8]. 

On  the  other  hand,  little  concern  has  been  shown  about  scaling  o  in  radar 
applications  [1],  if  the  model  is  a  metal  of  high  o,  since  it  is  felt  that  high  o 
models  have  only  a  small  surface  component  of  E,  and  differences  between  the 
surface  current  distribution  on  a  poor  conductor  and  that  on  a  perfectly  conducting 
target  have  been  believed  to  be  unimportant. 

However,  a  simple  argument  based  on  Equation  (19)  leads  to  quite  different 


predictions.  For  a  plane  wave  incident  on  a  physical  scale  model,  electric  and 
magnetic  energy  densities  are  equal  to  each  other  outside  the  conductor,  and  the 
radar  cross  section  (RCS)  comes  equally  from  both  electric  and  magnetic  scattering. 
Therefore,  if  the  ratio  of  scattered  H  fields  shrinks  with  respect  to  the  ratio  of 
E  fields  in  accordance  with  Equation  (19),  then  measured  values  of  the  RCS  of 
scaled  targets  should  approach  one-half  the  value  predicted  for  an  ideal  model. 

Actually,  magnetic  losses  are  twice  as  important  as  the  above  argument  would 
indicate.  Inside  th'^  conductor,  the  time  average  of  the  electrical  energy  density 
is  cE^/2,  while  that  associated  with  the  H  fields  is  oE*/2w,  so  their  ratio  is 
eu/o,  which  is  assumed  small  in  agreement  with  Equation  (26).  Therefore,  the 
energy  density  in  a  conductor  is  predominantly  magnetic,  so  good  conductors  are 
essentially  not  penetrated  by  E  fields  satisfying  the  condition  that  eio  <<  o. 

Thus,  instead  of  measured  values  of  the  RCS  of  incompletely  scaled  targets 
approaching  one-half  the  value  predicted  for  an  ideal  model,  the  diminished  surface 
currents  and  declining  H  fields  associated  with  a  a  not  increasing  fast  enough  with 
frequency  for  proper  scaling  lead  to  a  ratio  of  measured  values  of  the  RCS  to  that 
predicted  for  an  ideal  model  which  will  eventually  vanish  [26]. 

These  general  concepts  remain  valid  in  variable  electromagnetic  fields,  even 
if  dispersion  is  present  [27].  Because  of  the  continuity  of  the  tangential 
components  of  E  and  H,  the  normal  component  of  S  in  Equation  (8b)  is  continuous  at 
the  boundary  of  a  body,  and  the  interpretation  of  S  as  an  energy  flux  density 
carries  over  from  a  vacuum,  across  surfaces,  into  more  complicated  bodies. 

However,  in  the  presence  of  dispersion,  the  interpretation  of  energy  density  is 
more  complicated,  because  the  presence  of  dispersion  generally  signifies 
dissipation  of  energy,  since  a  dispersive  medium  is  also  an  absorbing  one. 
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SUMMARY  AND  CONCLUSIONS 

Since,  all  scattering  characteristics,  including  interference  and  diffraction, 
polarization,  and  creeping  and  traveling  (surface)  wave  phenomena,  should  be 
correctly  represented  by  rigorous  scale  modeling,  the  distribution  of  radiation 
scattered  from  a  full-size  target  is  expected  to  be  predictable  from  measurements 
on  a  faithful  scale  model. 

Therefore,  a  complete  and  rigorous  treatment  of  the  predictions  of 
electrodynamic  similitude  relating  scattering  from  scaled  and  full-size  targets  was 
developed,  starting  with  the  usual  assumption  that  e,  u,  and  o  do  not  change  with 
f.  Three  general  nonlinear  modeling  equations  in  six  variables  were  derived,  and  a 
complete  set  of  solutions  was  presented.  Scaling  laws  for  such  models  explicitly 
Involve  e,  m,  and  a.  However,  practical  difficulties  with  scaling  properties  of 


realistic  materials  now  exist  which,  for  most  scattering  measurements,  have 
discouraged  all  but  simple  geometric  scaling. 

Particular  emphasis  was  given  to  the  effects  of  differences  between  geometric 
and  complete  scaling  on  the  electromagnetic  fields  and  on  the  radar  cross  section 
(RCS),  and  the  effects  of  approximations  to  complete  scaling  were  evaluated. 
Conditions  were  obtained  on  properties  of  materials  required  for  models  made  from 
these  materials  to  accurately  simulate  systems.  Absorption  and  energy  balance  were 
also  treated,  and  the  influence  of  finite  o  on  surface  currents  was  shown. 


Scaling  laws  are  the  result  of  exact  conditions  on  each  coefficient  of 
Maxwell's  equations.  Since  the  claim  that  scaled-down  models  may  be  used  to 
represent  real  targets  is  based  on  linearity  of  Maxwell's  equations,  we  also 
examined  the  physics  involved  in  nonlocality  and  nonlinearity,  its  relation  to  the 
Maxwell  equations  in  macroscopic  media,  and  its  impact  on  physical  scale  modeling. 
Nonlocality  in  space  may  often  be  neglected  in  conductors  up  through  optical 
frequencies,  as  long  as  the  mean  free  path  for  conduction  electrons  is  small,  and 
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the  radar  wavelength  is  large,  compared  to  the  skin  depth.  When  these  conditions 
are  satisfied,  e,  m,  and  o  become  functions  only  of  f  in  dispersive  targets. 

The  behavior  of  electromagnetic  fields  was  also  examined  as  a  function  of  f. 

If  0  were  infinite,  as  is  usually  assumed,  satisfying  general  scaling  requirements 
would  be  elementary,  because  all  solutions  to  the  general  modeling  equations  are 
then  trivially  satisfied,  except  those  leading  to  geometric  scaling. 

Unfortunately,  it  is  physically  impossible  to  satisfy  scaling  conditions 
exactly  when  realistic  materials  are  used  and  the  change  in  f  is  appreciable, 
because  this  would  require  that  upper  bounds  on  physical  o  be  exceeded  at  higher  f, 
and  definite  limitations  also  exist  on  physical  g  and  p. 

This  paper  showed  that  incoruplete  scaling  of  apparently  nonabsorpti ve  targets 
may  have  profound  implications  for  model  measurements.  Both  electric  and  magnetic 
scattering  were  considered.  All  the  energy  density  in  the  reflected  wave  may  be 
lost  at  higher  f,  and  the  ratio  of  the  measured  RCS  of  geometrically  scaled  targets 
to  that  predicted  for  a  completely  scaled  model  will  eventually  approach  zero. 

RECOMMENDATIONS 

Even  though  we  have  concentrated  on  scattering  of  radar  waves  from  targets, 
very  few  conclusions  actually  depending  on  properties  of  waves  have  been  used, 
beyond  Equation  (20).  To  assess  effects  of  compromises  with  complete  scaling,  a 
more  realistic  calculation  of  the  precise  interaction  of  electromagnetic  waves  with 
physical  models  is  needed. 

First,  the  scattering  amplitude  should  be  calculated  for  electromagnetic  waves 
impinging  on  nondispersl ve  materials  In  the  proper  frequency  domain  with  a  more 
careful  description  of  fields  and  currents.  Such  a  description  is  quite  applicable 
up  through  microwave  frequencies  to  the  range  where  the  effects  of  electron  inertia 
begin  to  be  felt.  Both  reflection  and  absorption  coefficients  should  be 


calculated,  as  well  as  phase  angles.  Ratios  of  the  measured  radar  cross  section 


(RCS)  values  of  scaled  targets  to  those  predicted  for  an  ideal  nondisperslve  model 
should  also  be  calculated. 

These  calculations  should  be  extended  to  dispersive  materials  up  through 
higher  microwave  frequencies  into  the  range  where  the  effects  of  electron  inertia 
can  be  felt.  In  general,  as  f  increases  up  through  the  modeling  range,  and  becomes 
comparable  with  frequencies  corresponding  to  electron-lattice  interactions  and  to 
electron  motion  within  atoms  and  molecules  of  the  material  media  (optical 
frequencies),  absorption  takes  place,  which  leads  to  p,  e,  and  a  becoming  complex 
functions  of  f  with  real  and  imaginary  (reflective  and  absorptive)  parts. 

The  gradual  onset  of  absorption  in  good  conductors  as  f  increases  through  the 
modeling  range  permits  a  simple  description  in  terms  of  electron  inertia.  In 
conductors,  electron  inertia  is  usually  felt  at  a  lower  f  than  other  forms  of 
absorption,  and  is  also  the  first  cause  of  dispersion.  The  Drude  theory  describes 
these  dispersive  effects  of  electron  inertia  well  with  a  complex  e,  and  can  be 
applied  to  metals  as  o  decreases  with  increasing  f  through  the  modeling  range  from 
the  X-band  used  on  full-scale  ships  to  the  shortest  X  considered  for  physical  scale 
models.  Ratios  of  the  measured  RCS  values  of  scaled  targets  to  those  predicted  for 
such  a  dispersive  model  should  also  be  calculated  and  compared  with  those  obtained 
for  the  nondispersive  case. 
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